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We report experimental investigation of hyperfine optical pumping and near resonance light scat- 
tering from a high-density and ultracold atomic *^Rb gas. The atomic sample, having a peak density 
~ 5 • 10^'^ atoms/cm'^, temperature ~ 65 fjK and initially prepared in the F = 1 lower energy *^Rb 
hyperfine component, is optically pumped to the higher energy F = 2 hyperfine level. Measurements 
are made of the transient hyperfine pumping process and of the time evolution of scattering of near 
resonance probe radiation on the F = 2 —> F' — 3 transition. Features of the density, detuning, 
and probe laser intensity dependence of the signals are attributed to the high peak density and 
consequent large optical depth of the samples. 
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I. INTRODUCTION 

Development of techniques to cool and trap atomic 
gases [ll, Q] has revolutionized many traditional atomic 
physics research areas and, at the same time, has stim- 
ulated new connections and types of interdisciplinary 
specialization. Even if one does not consider vital re- 
search involving quantum degenerate gases 0], ar- 
eas such as quantum information 0, quantum optics 
0) B H) E ) precision measurements [111, [l^ . [TbII . plasma 
physics [3, , and molecular spectroscopy [l7[ have 
been transformed by the combination of ultracold experi- 
mental facility and theoretical understanding of the phys- 
ical processes. Among the important characteristics of 
utilized atomic gases are the gas temperature, density, 
and spin polarization. For example, storage of individual 
photon wave packets in ultracold gases combines quan- 
tum optical techniques of coherent dark state formation 
and electromagnetically induced transparency to address 
a critical area of quantum information processing [Tsl-fl^. 
The lifetime of the atomic spin wave, which determines 
the storage time of the photonic information in the form 
of a dark state polariton, depends on the gas temper- 
ature, collision rates, and local magnetic environment. 
In another research area, photoassociative formation of 
ultracold diatomic molecules depends quadratically on 
the density of the parent ultracold atomic gas. Likewise, 
by initiating formation of an ultracold plasma in a high 
density atomic gas, one can attain strong coupling with 
sufficient ionization of the gas. 

Over the past decade, a number of research groups 
[13, iH, m. 111 m, ii, ill have concentrated on near- 
resonance light scattering in ultracold atomic gases. A 
particular interest has been in coherent multiple scatter- 
ing of light when the influence of atomic motion is of 



minimal importance. Then the natural length scale for 
the photon multiple scattering is the optical mean-free- 
path, given by / = l/pa, where p is the atomic density 
and cr the cross-section for atomic light scattering in a 
weak field. For light of wave vector magnitude k = 2tt/X, 
a useful dimcnsionless parameter is kl [23|. When kl 
>> 1, the so called weak localization regime [1^, light 
scattering can be thought of as a sequence of scattering 
and propagation events 23]. In this regime, an observ- 
able that survives configuration averaging is the coherent 
backscattering cone, which is a few milliradian- width fea- 
ture that displays an enhancement of as much as a factor 
of two over the incoherent albedo for back scattered light 
[28 , 29, 30]. In ultracold atomic gases, the angular shape 
and peak enhancement of this spectral feature has been 
studied for a wide range of conditions, including spectral 
detuning from resonance, light polarization, probe light 
intensity [3l|, sample size and optical densitj^ and ex- 
ternal magnetic field [13, [H, H H 113, Hi S ES] ■ 
Observation of the coherent backscattering cone estab- 
lished that multiple light scattering, even for resonance 
radiation, is a coherent process. When the atomic den- 
sity is significantly increased, so that the parameter kl 
~ 1 (the so-called loffe-Regel condition [2?!), there are a 
number of atoms within a volume l/fc"^, and light scat- 
tering becomes a cooperative process [H [H [4^ . In 
this regime, it is expected that a number of fascinating 
quantum optical processes, including Anderson localiza- 
tion of li ght I44l.[45ll46ll47l . l48l |. and random atom-based 



lasing [49|, |50|, |51 
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, 52I |53||. mav emerge. It is these pro- 
cesses, and the conditions under which they might be 
experimentally studied, that are the main motivation for 
the current research program. 

Before we present our experimental procedures, results 
and associated discussions, let us make a few remarks 
concerning the definitions and parametrization we fol- 
low throughout the paper. Most of the reported exper- 
imental data relate to a relatively high density atomic 
system. For the considered densities, which are at the 
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level of more than one atom per cubic wavelength, it 
appears that the atomic system should have properties 
approaching those of a bulk medium. From a classical 
point of view the relevant theoretical description then 
would require us to introduce a mesoscopically smoothed 
dielectric constant e and treat the light scattering pro- 
cess through secondary or multiply scattered waves cre- 
ated by fluctuations in e. If the fluctuation scale of e 
were much smaller than its mean value then the light 
transport through the sample would be fairly approxi- 
mated by a perturbation theory analysis of the macro- 
scopic Maxwell equations, [sj]. Apparently, such an ap- 
proximation would be excellent in a weak disorder regime 
[23l | at low densities when e ~ 1. For the peak values of 
the atom densities attained in our experiments the vari- 
ation of the dielectric constant is expected to be still 
quite close to the vacuum level [l^. This justifies the 
parametrization of our data in the standard terms nor- 
mally valid for a dilute configuration, such as the vacuum 
scattering cross section, the mean free path for light scat- 
tering, and the optical depth. These quantities seem to 
us to be convenient for qualitative description and dis- 
cussion. However at certain points of our analysis we will 
pay attention to the evidence of the effects of macroscopic 
electrodynamics, see Section III. C. We would also like to 
point out that at a further increase in density by only sev- 
eral times, the situation dramatically changes and turns 
toward a highly disorder regime. As expected, in this 
case the average dielectric constant would have the same 
order of magnitude as its fluctuations, and the system 
would demonstrate some specific spectral properties of a 
collectivized quantum system ji^lisj. These properties 
have a signature of strong correlation with each particu- 
lar configuration of atomic scatterers, and can be linked 
with the Anderson-type localization conditions for elec- 
tromagnetic waves. Part of the single dipole excitations 
associated with the collective resolvent poles would have 
a long-lived sub-radiant nature. 

In the present paper, we describe our experimental ap- 
proach to obtaining atomic conditions, and particularly 
atomic densities, in the regime where k/ ~ 1. We also de- 
scribe the method we use to prepare the atomic sample so 
that large orders of multiple scattering may be obtained, 
along with the associated experimental observations and 
their interpretation. This is followed by presentation and 
discussion of experimental results of density, detuning, 
and light intensity dependence of light scattering on the 
nearly-closed F = 2 — > F' 3 hyperfine component of 
the *^'^Rb D2 transition. 



II. EXPERIMENTAL APPROACH 

A schematic of the experimental apparatus is shown in 
Fig. 1, and the main optical excitation transitions used 
in the present experiments are shown in Fig. 2. In Fig. 1, 
the main focus of the instrumentation is an optical dipole 
trap formed in the focal region of a carbon-dioxide (CO2) 
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FIG. 1: Schematic drawing of the experimental apparatus. In 
the figure QUEST stands for quasi electrostatic trap, PMT 
refers to a photomultiplier tube, and CCD represents a charge 
coupled device camera. MOT refers to a magneto-optical 
trap, while CO2 laser indicates a carbon-dioxide laser. Draw- 
ing not to scale. 




FIG. 2: Schematic energy level diagram for Rb showing the 
optical transitions important for the experiments reported in 
this paper. Not to scale. 



laser. A 40 MHz acousto optical modulator (AOM) is 
used to switch and to direct approximately 50 W of the 
100 W CO2 laser output to the trapping region. The 
trapping beam is focussed to a radial spot size of ~ 55 
fim, and associated Rayleigh range of z^f 750 /im. Be- 
cause the operating trap wavelength of 10.6 fj,m is far 
longer than those of the Rb resonance transitions, the 
dipole force is closely proportional to the static dipole 
polarizability of atomic Rb; the trap is then referred to 
as a quasistatic electric dipole trap (QUEST) 0. The 
trap depth is approximately 650 ^K, and has associated 
measured trap angular frequencies lo^ = 2tt ■ 50 rad/s 
and = 27r • 1300 rad/s. As described in detail else- 
where [5^, the QUEST is loaded with atoms that have 
been collected from a surrounding thermal Rb vapor and 
cooled in an overlapping magneto optical trap (MOT). 
The QUEST atoms are initially loaded into the lower en- 
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ergy F = 1 hyperfine component of the Rb ground level. 
Following QUEST loading, the MOT lasers and magnetic 
fields are shut off, and the atoms in the dipole trap colli- 
sionally thermalize to a temperature of approximately 65 
fiK. After this natural thermalization process, approxi- 
mately 15 % of the atoms originally in the MOT have 
been transferred to the QUEST. Measurements of the 
QUEST characteristics, after the hold period, by absorp- 
tion imaging, parametric resonance, and the measured 
number of atoms transferred show a sample with peak 
density about 5 ■ 10"'^^ atoms / cm^ and a temperature of 
To — 65 ijlK. The 1/e lifetime of the confined atoms is 
greater than 5 s, limited by background gas collisions. 

The primary goal of the experiments is to study light 
scattering on the nearly closed F — 2 ^ F' — ?> hyperfine 
transition associated with the D2 resonance line. How- 
ever, the atoms are initially loaded into the QUEST in 
the F = 1 lower energy hyperfine component, and must be 
optically pumped into the higher energy F = 2 level. This 
is accomplished by the MOT repumping laser beams, 
tuned resonantly to the F = 1 — > F' = 2 transition 
as indicated in Fig. 2. The repumper laser is an external 
cavity diode laser, and is locked to an *^Rb saturated 
absorption feature associated with the F — \ F' — 2 
hyperfine transition. The laser bandwidth is approxi- 
mately 0.5 MHz. The repumper delivers a beam of max- 
imum intensity ~ 4 mW/crn? and is directed along the 
same optical paths as the trapping laser beams. Re- 
pumper switching is controlled with an acousto optical 
modulator. There are three counter propagating pairs 
of such beams directed towards the sample along three 
orthogonal directions. The repumper intensity of 4 
mW/ cm? corresponds to an on- resonance saturation pa- 
rameter larger than unity. With reference to Fig. 1, the 
light scattered out of the repumper beams by the atom 
sample is collected at an angle of 45 degrees from the hor- 
izontal pair of repumper beams and at 90 degrees from 
the vertical beam pair (not shown). Approximately 1 % 
of the resulting scattered light is collected with a field 
lens and launched into a multimode fiber which, in turn, 
transports the light to an infrared sensitive and refriger- 
ated photomultipher tube (PMT) operating in a single 
photon counting mode. The photon counting pulses are 
amplified with a fast amplifier and then binned in a mul- 
tichannel scalar (MCS) with a time resolution of 5 ns. 
Depending on the counting rate, a data record of several 
thousand experimental realizations is necessary to obtain 
sufficient counting statistics. 

In the main experimental protocol, a probe beam 
tuned in the spectral vicinity of the F = 2 — > _F' = 3 
nearly closed transition is directed towards the sample, 
and the resulting scattered light signals collected by the 
same optical arrangement described above. The probe 
laser is of the same design as the repumper laser, has a 
bandwidth ^ 3 MHz, and is switched and directed by 
an acousto optical modulator towards the sample. Be- 
cause of constraints on the vacuum chamber geometry, 
the linearly polarized probe beam is directed (see Fig. 
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TABLE I: QUEST parameters relating the peak transverse 
optical depth on the F — 1 F' — 2 transition to the peak 
sample density and the Gaussian radii of the atomic cloud. 

1) at an angle of approximately 30 degrees away from 
the fluorescence collection direction. The probe beam 
is also directed downwards at an angle of 30 degrees. 
The collection and electronic accumulation of scattered 
light signals is the same as with the repumper signals. 
We should point out that the two protocols described 
above are applied separately or together, and with var- 
ied relative time separation, depending on the specific 
experiments of interest. These different approaches are 
described individually, and as appropriate for each data 
set, in the following sections. 

To close this section, we point out that in some of 
the experiments reported here the atomic density was 
varied over a wide range. This was accomplished by al- 
lowing for a period of ballistic expansion of the cloud 
after the QUEST was turned off. The atomic sample 
temperature is known, so this procedure allows the peak 
density or the peak optical depth to be determined. As 
the sample is well approximated by a two-axis Gaussian 
atom distribution '50|, the two Gaussian radii and the 
peak atom density (or the total number of atoms in the 
sample), are sufficient to determine the two peak opti- 
cal depths characterizing the sample. We summarize in 
Table 1 the peak transverse optical depth bt the peak 
atom density at the center of the sample Uo, the trans- 
verse Gaussian radius Tq, and the longitudinal Gaussian 
radius Zo- The optical depth refers here to that of the 
nearly closed F — 2 ^ F' — 3 hyperfine transition, and 
is obtained by using a resonance scattering cross section 
of 1.36 X 10~^cm^. The optical depth dependence may 
readily be rescaled to other F — > F' transitions [s^ . 

III. RESULTS 

A. Hyperfine Optical Pumping 

In this section, we focus on the hyperfine optical pump- 
ing process, by which the ®^Rb atoms in the QUEST are 
brought from the F = 1 to the F = 2 ground hyperfine 
component. Fluorescence signals taken over a three order 
of magnitude range of sample density are shown in Fig. 
3. These densities, as indicated in the figure legend, rep- 
resent the peak density in the central region of the atomic 
sample. As described in an earlier section of the paper, 
the sample density is reduced from its maximum value 
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FIG. 3: Time-dependence of the light scattered during the 
F = 1 F' = 2 optical pumping process. From upper to 
lower, the curves correspond to decreasing atomic density; the 
values refer to the peak density at the center of the sample. 
The correspondence between the peak density, the transverse 
optical depth, or the sample dimensions can be found in Table 
1. 



by allowing for a selected period of free expansion of the 
sample prior to switching on the repumper laser and col- 
lecting the resulting scattered light signals. Further ex- 
pansion of the sample during the data acquisition period 
is negligible. Several features of the data are apparent 
from Fig. 3. First, there is a small background signal of 
a few counts remaining even after the transient signals 
have decayed away. This behavior is more apparent in 
the lower density data, which has an associated faster de- 
cay rate than the higher density data. This signal is due 
to hot atom fluorescence excited by the intense repumper 
laser beams. Second, the decay of the signals is well ap- 
proximated by a single exponential time decay, with the 
higher density data decaying at a much slower rate than 
the lower density measurements display. Finally, the to- 
tal integrated intensity does not seem to vary much as 
the density is varied. These three features are summa- 
rized in Fig. 4. In this figure, the measured quantities 
are shown as a function of the peak transverse optical 
depth bt- The connection between the various values of 
bt and the peak atom density can be found in Table 1. 
In Fig. 4, we see that the peak scattered light intensity 
and the total (time-integrated) intensity does not vary 
strongly with optical depth. However, the exponential 
decay rate associated with the pumping process depends 
approximately exponentially itself on bt- The time scale 
associated with the optical pumping process may be es- 
timated by considering that the repumper laser beam is 
intense enough that it saturates the repumper transition. 
Then the time scale for the pumping process transferring 
an atom from the lower (F=l) to the upper (F = 2) hy- 
perfine component is approximately 2to, where Tq ^ 26 ns 
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FIG. 4: Variations, with peak transverse optical depth, of 
measured quantities associated with the time and density de- 
pendence of the optical pumping transients, as extracted from 
the data of Figure 3. (a) Integrated and peak scattering sig- 
nals, and (b) decay rate of the repumper transients. 



is the radiative lifetime of the excited state. On the aver- 
age, then the repumper beam penetrates into the sample 
a distance of one optical mean free path ^ = ^ in this 
time. Here p is the average atom density and a is the 
optical scattering cross section For a sample of size 
To, an estimate of the pumping time is then T ^ 2to^. 
This simple formula provides estimates in fair agreement 
with the results of Fig. 4(b). 

Finally, we point out that in spite of the apparently 
simple phenomenology of the optical pumping dynamics, 
the physical process is in reality quite complex. For ex- 
ample, the atom sample is initially in the lower energy F 
= 1 hyperfine level, and after the process is terminated, 
is nearly completely pumped into the F = 2 level. At 
the same time, the atom sample is initially very optically 
deep to the repumper laser and optically thin to the F = 2 
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^ F' = 2 inelastic Raman decay. These roles are reversed 
but in a spatially inhomogeneous way as the pumping dy- 
namics take place. Beyond this, multiply scattered light 
on both main transitions should participate significantly 
in the entire dynamics (this is ignored in the estimate of 
the pumping time above). Finally, the entire process is 
made yet more complex by the inelastic components in 
the scattered light, these being generated by the intense 
repumper laser. As this process is not the main focus 
of the present work, we defer theoretical modelling and 
more comprehensive experiments to a later study. 



B. Combined Hyperfine Optical Pumping and 
Light Scattering on the F = 2 ^ F' = 3 Transition 
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FIG. 6: Steady state signals as a function of F = 2 ^ F' = 
probe laser intensity, as extracted from the data of Fig. 5. 
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FIG. 5: Probe intensity dependence of the combined tran- 
sients due to hyperfine optical pumping on the transition and 
light scattering on the nearly-closed F = 2 ^ F' = 3 transi- 
tion. 

We now consider measurements made of the time evo- 
lution of the light scattered from the sample when the 
repumping laser and a second probe laser beam are ap- 
plied simultaneously to the sample. One main goal of 
these measurements is to find the approximate intensity 
range suitable for optical probing on the F = 2 ^ F' = 3 
transition, and to discover if there is significant trapping 
of probe light within the sample through the combined 
action of the repumper and probe excitations. For these 
measurements, the probe laser is tuned resonantly to the 
nearly-closed F ^ 2 ^ F' ^ 3 hyperfine transition as- 
sociated with the D2 component of the 5s - 5p multi- 
plet. Results from these experiments are shown in Fig. 5, 
where the different curves correspond to different probe 
laser intensities. The repumper intensity is fixed, and 
at the same level as the results shown in Fig. 3. We see 
from Fig. 5 that the time evolution of the signals depends 
strongly on the probe laser intensity. For lower intensi- 



ties, the time evolution is similar to that of the repumper 
laser alone, while for greater probe intensities, the signals 
more rapidly reach a nearly steady state level. The re- 
sults may be qualitatively understood by considering that 
the atomic sample is initially optically very deep on the 
repumper transition, but becomes increasingly less so as 
a function of time. This behavior is evident from the data 
of Fig. 3. On the other hand, as the atoms are pumped 
to the F = 2 level, the optical depth associated with the 
F = 2 ^ F' = 3 probe transition becomes increasingly 
large, the transverse optical depth becoming similar to 
that of Table 1 for later times. As the probe transition 
is nearly optically closed, we expect a steady scattering 
rate to be reached, and the data of Fig. 5 supports that 
expectation. The steady state level depends on the probe 
laser intensity, and the variations of this level are summa- 
rized in Fig. 6. There we see that the steady state level 
increases approximately linearly for lower intensities, but 
saturates at higher levels. However, the saturation pa- 
rameter of unity for the probe transition is reached for 
a probe intensity of Is ^ l.lmW / cm?; the effective sat- 
uration evident in Fig. 6 occurs for much higher probe 
intensities than that. This may be physically understood 
through the fact that most of the atoms contributing to 
the signal are deep inside the atomic cloud, where the 
probe intensity is sharply decreased by scattering out of 
the coherent beam to be much smaller than the peak 
external value. 
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C. Probe Light Scattering on the F = 2 ^ F' = 3 
Transition 

1. Density Dependence 

We now turn our attention to the main focus of this 
report, which is the situation where the hyperfine op- 
tical pumping process is completed prior to initiation 
of probe light scattering in the spectral vicinity of the 
F = 2 — > F' = 3 transition. In this case, the function of 
the hyperfine optical pumping is to transfer the atoms to 
the F = 2 ground level hyperfine component; when this 
process is complete, the peak resonance optical depth on 
the probe transition is about 165, as indicated in Table 
1. We first consider the resonance response as a function 
of peak atom density under weak-field probe conditions 
of ^ 630 jiW/cm? (see Fig. 6). For these measurements, 
the atom sample is exposed to a nearly rectangular tem- 
poral pulse of 2 fjLS duration; this pulse has a 20 dB rise 
and fall time of about 100 ns. As shown in Fig. 7(a) and 
7(b), the scattered light transients for all densities consist 
of a rather rapid increase to a steady level, followed by 
a several hundred nanosecond temporal decay after the 
probe laser beam is extinguished. 

We further analyze these results by considering the in- 
tegrated probe signal (as in Fig. 7) as a function of de- 
creasing atomic density. However, instead of the atomic 
density, we parametrize the dependence in terms of the 
peak transverse optical depth bt through the center of the 
ellipsoidal atomic sample. The optical depth is the nat- 
ural parameter to describe many characteristics of light 
scattering and diffusion in dense scattering media, and 
the results described here are no exception. Note that 
there are nominally two ogtical depths required to de- 
scribe our atomic sample [56| . These are the peak trans- 
verse optical depth bt, as just mentioned, and the lon- 
gitudinal optical depth which is typically more than 
10 times larger than bt. The variation of the integrated 
scattered light signals with bt is shown in Fig. 8, where 
it is seen that the signals increase as the optical depth 
decreases. This behavior can be physically understood 
from the fact that for very large bt, only the atoms near 
the sample surface contribute to the scattering signal. 
However, as the density, and correspondingly the optical 
depth, decreases, more of the atoms in the sample partic- 
ipate in the scattering, and the resulting signal increases. 

To more quantitatively understand these results, we 
consider a simple model of light scattering from a spa- 
tially Gaussian atom distribution that is entirely illumi- 
nated by the probe beam. This model shows that the 
total (spatially integrated) scattered light intensity Ps 
should vary with optical depth b as [5^ 
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FIG. 7: Variations of the F = 2 ^ F' = 3 resonant probe 
response with density and time. Note that the maximum 
intensity decreases with increasing density, at constant atom 
number. Note also the emergence of a long-time decay after 
the probe excitation is shut off. (a) Dependence for higher 
densities, (b) Dependence for lower intensities. 



the incident probe intensity, assumed to be uniform over 
the cross-section of the sample, N is the number of atoms 
in the sample, and ctq is the resonance scattering cross 
section, b is the peak optical depth. 

Here we make the assumptions that the optical depth 
corresponds to bt (viz. Table 1) and that 77 is an empirical 
scaling parameter for the optical depth that is meant to 
measure the actual nonuniform spatial density distribu- 
tion of the atomic sample and the limited penetration of 
the probe beam into the sample volume. This heuristic 
model describes the data of Fig. 8 well with values of A 
= 25(2) Kcounts, and 77 = 0.08(2). The small value of 77 
is consistent with our expectation that the probe beam 
does not significantly penetrate the atomic vapor, even 
for the smallest transverse optical depth of about 5. The 
result of the two parameter fit of the expression for Ps{b) 
to the density dependent data is shown as a solid curve 
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in Fig. 8. 

Although this model gives a decent fit to the exper- 
imental data, it does not take into account the double 
Gaussian spatial distribution of atoms in the cloud or 
the angular distribution of the scattered intensity. More 
importantly, for the higher atomic densities correspond- 
ing to the larger optical depth parts of the data set, the 
scattering cross section departs from the free atom ex- 
pression. As we describe in detail elsewhere , we have 
performed theoretical analysis and Monte Carlo simula- 
tions of the data in Fig. 8, and of data presented later in 
this paper. The fully quantal approach used has been de- 
scribed in detail previously , but has been generalized 
to account for the bi-Gaussian sample shape and for the 
particulars of the experimental excitation and detection 
geometries. 
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FIG. 8: Variations of the integrated F — 2 ^ F' = 3 reso- 
nant probe response with peak transverse optical depth. Note 
that the maximum intensity decreases with increasing optical 
depth, at constant atom number. 
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FIG. 9: Representative spectral variations of the transient 
light scattering response associated with the F = 2 — > F' = 
3 probe transition. Positive (higher frequency) detunings are 
shown in (a), while negative (lower frequency) detunings are 
shown in (b). These data are recorded under conditions of 
the peak transverse optical depth bt = 165. 



2. Detuning Dependence 

We now turn to the detuning dependence of the light 
scattering signals. In these experiments we have mea- 
sured the temporal response of the sample to a, 2 fis 
probe pulse for detunings A in the range ± 24 MHz 
from the atomic F — 2 ^ F' = 3 resonance. Positive 
values of A correspond to probe frequencies larger than 
the isolated atom resonance frequency. Representative 
measurements are given in Fig. 9 (a), (b), where we see 
that the response to the nearly rectangular probe pulse 
consists of a quite rapid increase in signal to a steady 
value (within the measurement statistics), followed by a 
clear decaying signal extending for nearly a /is after the 
probe beam is extinguished. These data are recorded un- 
der conditions of the peak transverse optical depth bt = 



165. In analysis of this data, we first consider the tem- 
porally integrated signals as a function of detuning from 
atomic resonance; the resulting spectral variations are 
shown in Fig. 10. There we see that the response is ap- 
proximately spectrally symmetric and has a full-width at 
half-maximum of about 24 MHz, corresponding to four 
times the natural width of the F — 2 F' — 3 tran- 
sition under weak probe field conditions. Qualitatively, 
we expect broadening of the atomic response due to the 
large optical depth of the sample. We also note a small 
shift of the resonance response to frequencies lower than 
the single atom resonance. The frequency shift of about 
-0.8 MHz is of the right magnitude to correspond to the 
well-known Lorentz-Lorenz (local field) shift [5^. How- 
ever, further measurements at higher density and with a 
spectrally narrower probe would be necessary to quanti- 
tatively examine that possibility. 
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FIG. 10: Spectral variation of the temporally-integrated scat- 
tered light intensity on the F = 2 ^ F' = 3 transition. In 
this data the repumping process occurs first, and the near- 
resonance scattering signals are recorded after the repumping 
laser is turned off. These data are recorded under conditions 
of the peak transverse optical depth bt — 165. 



To analyze further this data, we return to the heuristic 
model we used m Section III.C.l, but taking into account 
the fact that the sample size remains fixed in these mea- 
surements, and that the optical depth varies with spectral 
detuning from atomic resonance. This model then gives 



Ps{b) = Aj2{-iy 



-1 ivbr 



(3.2) 



n=l 



value of the empirical parameter rj is also in satisfactory 
concordance with that obtained from fitting the density 
dependent data to the same model. The laser spectral 
width used in the fit is 3 MHz, consistent with the mea- 
sured probe spectral profile. 

Finally, we consider the longer time transients that ap- 
pear once the probe laser beam is switched off. In the 
present case, when the probe is extinguished on a time 
scale fast compared to the emergence of light from the 
sample, the time scale is determined by diffusion of light 
from the interior of the sample to the surface, where it 
escapes and may be detected. We present in Fig. 11 
an extended experimental run of longer than 10 hours. 
In Fig. 11(a) the entire response, which consists of the 
hyperfine pumping process and the probe excitation is 
displayed. In this case the probe pulse is 5 fis. In the fig- 
ure, we see that there is a significant long time transient 
after the probe is shut off (vertical line). This long-time 
region is displayed with expanded horizontal and verti- 
cal scales in Fig. 11(b), where it is seen that the decay is 
characterized by multiple time scales, and that the decay- 
ing signals persist for longer than 2 ^s. A fit of this data 
to a double exponential decay resulted in a high quality 
representation of the data with time constants of 100(10) 
ns and 500(50) ns. This fit is shown as the dashed fine 
in Fig. 11(b). 

To put the longer of these time scales in context, we 
again consider a spherically symmetric Gaussian ultra- 
cold atomic cloud. For that case, and for large optical 
depths, the longest lived diffusive mode has a decay time 

m 



where 



36^ 
5.35^2^° 



(3.4) 



Here A is proportional to lofg, where Iq is the probe 
intensity, assumed uniform over the cloud dimensions and 
To is the Gaussian radius of the cloud. The quantity bo 
is the peak optical depth, while S accounts for the ob- 
served shift of the response profile. The spectral width 
7' is assumed to be the 6 MHz width of the isolated atom 
response plus the spectral width of the probe laser. This 
part of the model assumes that the spectral response is 
the convolution of two Lorentzian line shapes. As before, 
we take bo to be the peak on-resonance optical depth 
through the center of the cloud, and use 77 and d as fit- 
ting parameters. The result of a fit with values of 77 = 
0.10(1) and S = 0.8 MHz is shown as a solid curve in 
Fig. 10. In this fit, the overall scattered intensity is 
scaled by the parameter A, which measures the fraction 
of the total intensity that is collected and measured in 
the experiment. The result is quite satisfactory, and the 



where b is the peak optical depth and Tq is the ^ 26 
ns natural decay time of the excited level. If the atomic 
cloud in the present experiments were uniformly optically 
excited, the longest decay time should be about 40 /is, 
much longer than what is observed. We attribute this pri- 
marily to the relatively limited penetration of the probe 
beam into the atomic sample under the conditions of the 
present experiments. Using the same model considered 
before, if we replace b with rjbt = 16.5, this formula pre- 
dicts a decay time of 420 ns, consistent with what we 
have measured. 

We should point out that other factors, including sup- 
pression of the atomic cross section and the important 
role of atomic motion is also not negligible in these es- 
timates. Motion of the atomic scatterers, for example, 
leads to a random walk in frequency space for the mul- 
tiply scattered photons and, even for the atomic sample 
temperature of 65 , leads to a shortened lifetime for 
the longest lived diffusive mode. 
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FIG. 11: Transient light scattering response associated with 
the F = 2 ^ F' = 3 transition. In this data the repumping 
process occurs first, and the near-resonance scattering signals 
are recorded after the repumping laser is turned off. Both the 
repumping signal and the overall transient response is shown 
in part (a). The vertical line represents where the probe beam 
is turned off. In part (b) we expand the scales to show the 
long-time transient after the probe beam is turned off. 



IV. CONCLUSIONS 

In this paper we have reported studies of near reso- 
nance hght scattering from dense and ultracold atomic 
rubidium. The density is sufficiently large that the ex- 
perimental conditions are close to satisfying the loffe- 
Regal criterion for localization on the optically closed 
F = 2 ^ F' — 3 transition and on the optically open 
F = I ^ F' = 2 hyperfine repumper transition. How- 
ever, for the _F = 2 — > F' = 3 transition, external the 
large optical depth of the sample does not permit injec- 
tion of optical excitation to the highest density central 
portion of the sample. Observation of localization ef- 
fects then do not seem possible using this scheme. For 
the F = 1 ^ F' = 2 repumper transition, there is an 
open F' — 2 ^ F' — 2 inelastic Raman decay channel 
which strongly inhibits multiple scattering, again making 
it challenging to observe localization effects. In spite of 
this, optical excitation on these transitions has revealed 
a variety of optical pumping and dynamical propaga- 
tion effects which depend significantly on detuning from 
atomic resonance and on the optical depth of the atomic 
medium. For the case of a resonant probe beam incident 
on the cloud at maximum density, the light minimally 
penetrates into the medium, and the observed effects are 
due to scattering from atoms within a few optical mean 
free paths of the sample surface. This observation is sup- 
ported by the relatively short diffusive lifetime of a few 
hundred ns for light emerging from the sample after the 
probe beam is turned off. Finally, to insure that light 
emerges from the sample via the longest lived mode will 
require alternative experimental schemes whereby pho- 
tonic excitation may be injected directly into the highest 
density portions of the sample. 
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